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ABSTRACT 


Paramterizations of single nucleon removal from the electromagnetic and strong 
interactions of cosmic rays with nuclei are presented. These parametrizations are based 
upon the most accurate theoretical calculations available to date. They should be very 
suitable for use in cosmic ray propogation through interstellar space, the Earth’s 
atmosphere, lunar samples, meteorites, spacecraft walls and lunar and martian habitats. 


PACS: 97.70.S, 96.40.De 



I. INTRODUCTION 


Galactic cosmic rays (Shapiro 1983, Friedlander 1989) are very high energy 
particles confined to the region of our Milky Way galaxy. They consist of about 98% fully 
stripped nuclei including protons and alpha particles, and about 2% electrons and positrons 
(Simpson 1983). Of the nuclear component, about 87% is hydrogen, about 12% is helium 
and the other 1% consists of heavier nuclei. Of these heavier nuclei, the CNO group and Fe 
are the most abundant with a typical energy of about 1 GeV/N. Even though these heavy 
nuclei are not very abundant, they are very penetrating due to their large mass and high 
speed. 


An understanding of the interactions of galactic cosmic ray nuclei with not only 
hydrogen and helium but also with heavier nuclei is important for several reasons: 

!. Knowledge of the cosmic ray spectrum at the top of the Earth's atmosphere and 
knowledge of the composition of the interstellar medium and heliosphere enables one to 
determine the cosmic ray spectrum at the source (Simpson 1983). The interstellar medium 
(Field 1986) consists primarily of hydrogen and helium so that cosmic ray interactions with 
these nuclei are the most important (Austin 1981; Ferrando et al 1988). However carbon, 
nitrogen and oxygen are also present in the interstellar medium (Morton 1975; Karttunen, 
Kroger, Oja, Poutanen and Donner 1987) and one anticipates that the understanding of 
cosmic ray interactions with these heavier nuclei may be needed in the future. 

2. Knowledge of the spectrum at high altitude and knowledge of the composition of the 
Earth's atmosphere enables one to determine the cosmic ray spectrum at the top of the 


atmosphere (Wilson, Townsend and Badavi 1987a). 


3. The radiation environment inside a spacecraft, due to solar and galactic cosmic rays may 
be determined (Wilson and Townsend 1988). Such knowledge is important for lunar and 
martian habitats and other long duration space flights (National Council on Radiation 
Protection and Measurements 1989; Joselyn and Whipple 1990; Rester and Trombka 
1989). 

4. Studies of the history of extraterrestrial matter (such as lunar samples, meteorites and 
cosmic spherules and dust found in deep sea sediments) and also of the history of cosmic 
rays themselves can be made with the knowledge of the production rate of various nuclides 
(Reedy 1987; Reedy, Arnold and Lai 1983). 

5. Cross section parametrizations of cosmic ray nuclei interacting with arbitrary target 
nuclei including those targets heavier than helium are required in the interpretation of 
emulsion data and in the the interactions of cosmic rays with air. (Gaisser 1990; Gaisser, 
Stanev, Freier and Waddington 1982; Gaisser and Stanev 1983, Shapiro and Silberberg 
1970). 


The basic nucleus-nucleus interaction that a cosmic ray undergoes can occur 
mainly via the Strong or Electromagnetic (EM) force. (Actually the study of nucleus- 
nucleus collisions began in cosmic ray studies (Goldhaber and Heckman 1978; Bradt and 
Peters 1948, 1949, 1950; Kaplan, Peters, Reynolds, and Ritson 1952).) Strong interaction 
processes (Goldhaber and Heckman 1978; Gyulassy 1981; Benesh, Cook and Vary 1989) 
have been studied extensively and quite recently the study of Electromagnetic processes in 


high energy collisions has begun (Bertulani and Baur 1988). 


To study the propagation of cosmic rays through interstellar space, the Earth s 
atmosphere or a spacecraft wall it is not enough to have a good understanding of the 
nucleus-nucleus interaction cross section as input to a transport computer code (Wilson, 
Townsend, Schimmerling, Khandelwal, Khan, Nealy, Cucinotta, Simonsen, Shinn and 
Norbury 1991). These codes can be very complex and therefore require simple expressions 
for the cross sections rather than the use of large data bases or complicated theoretical 
models (Wilson and Townsend 1988). Thus there has been a considerable effort to 
parameterize the cross section expressions so that the only required inputs are the nuclear 
energies and charge and mass numbers (Letaw, Silberberg and Tsao 1983; Silberberg and 
Tsao 1973, 1990; Townsend and Wilson 1986; Norbury, Cucinotta, Townsend and 
Badavi 1988; Wilson, Townsend and Badavi 1987a,b). 

In order to understand cosmic ray transport through the interstellar medium, the 
early work on parametrizations (Rudstam 1966; Letaw, Silberberg and Tsao 1983; 
Silberberg and Tsao 1973, 1990) concentrated primarily on proton-nucleus interactions due 
to the fact that the interstellar medium consists primarily of hydrogen. However based on 
the 5 items listed above it would also be very useful for a wide variety of cosmic ray 
studies to have accurate parametrizations for any nucleus-nucleus interaction. It is the aim 
of the present work to provide such a parametrization. Actually such parametrizations 
(Wilson, Townsend and Badavi 1987b) have already been formulated and give good 
results for the removal of many nucleons. However for removal of only a few nucleons 
from heavy nuclei, the parametrizations (Wilson, Townsend and Badavi 1987b) sometimes 
give poor results. In fact a whole new approach to the parametrization of few-nucleon 
removal cross sections in nucleus-nucleus interactions is required. In the present paper an 


accurate parametrization of single-nucleon removal cross sections is presented. Future 
work will discuss the removal of more nucleons. When this program is completed we will 
have available accurate parametrizations of few-nucleon removal cross sections in nucleus- 
nucleus interactions. When combined with the many-nucleon removal paramterizations 
(Silberberg, Tsao and Shapiro 1976; Wilson, Townsend and Badavi 1987b) and proton- 
nucleus parametrizations (Letaw, Silberberg and Tsao 1983; Silberberg and Tsao 1973, 
1990), there will be available accurate cross section parametrizations for arbitrary cosmic 
ray species interacting with arbitrary media. See also the work of Webber et al (1990). 

One approach to the parametrization of cross sections is to simply take all the 
available experimental data and fit a curve through it. However such an approach often 
requires a large number of adjustable parameters and may not be applicable to regimes 
where experiments have not been performed. A much more satisfying approach is to base 
one's parametrization on a physical theory or model that successfully describes the 
experimental data as well. This will be the approach of the present work. The various 
models and theories that have been developed will be collected together and parameterized. 
The whole method will require only one adjustable parameter (Xd in equation 29). 
Furthermore this parameter is not essential. Good results are obtained without it. It is only 
introduced to provide some fine tuning. 

A preliminary parametrization of the EM process has already been presented 
(Norbury, Cucinotta, Badavi 1988), which utilizes the Weizsacker-Williams (WW) method 
of virtual quanta (Bertulani and Baur 1988; Jackson 1975). However, since then the theory 
has been improved to include the effects of both electric dipole (El) and electric quadrupole 
(E2) interactions (Bertulani and Baur 1988; Norbury 1990a,b), which will henceforth be 
referred to as multipole theory in contrast to WW theory. In addition Benesh, Cook and 


Vary (1989) have recently provided a parametrization of the strong interaction single 
nucleon removal cross section. 


IT. STRONG INTERAC TION PARAMETRIZATIQM 

The parametrization due to Benesh, Cook and Vary (1989) is 

o(N) = x °G p esc < la > 

A 

for single neutron removal where N is the number of neutrons and A is the number of 
nucleons and 


o(Z) = Ja G Peso < lb > 

for single proton removal where Z is the number of protons. See also Norbury and 
Townsend (1990). Cg is the reaction cross section given by 

o G = 2 tt (b c - 4J 1 ) Ab (2) 

where 

Ab = 0.5 fm (3) 


and the critical impact parameter for single nucleon removal is 


b c = 1.34 fm [Ap 3 + Aj 73 - 0.75(Aj> 173 + Aj 173 )] 


(4) 


with Ap and At being the projectile and target nucleon numbers respectively. The single 
nucleon escape probability is 


Pesc = (l-f) + fexp-v 


(5) 


with 

f = 2 (1 " ^OS ®max) 


( 6 ) 


and 


sin 0max — 


be - Ab 
be 


(7) 


and 

v = Aq NN (8) 

7tb c 2 

where A is the nucleon number of the nucleus from which the nucleon is being removed, 
and be is the critical impact parameter for the single nucleon escaping and is given by b c in 
equation (4) but with Aj = 1 (or Ap = 1) if the nucleon is escaping from the projectile (or 
target). Thus the escape probability is independent of Aj (or Ap) as one would expect. CTnn 
is the nucleon-nucleon cross section which has been parameterized as (Wilson, Townsend, 
Nealy, Chun, Hong, Buck, Lamkin, Ganapol, Khan and Cucinotta 1989) 

Onn = (1 +=^— ! ) {40 + 109 cos(0.199tiVE/180) exp[-0.451(Ti a b - 25) 0 258 ] } mb (9) 

Mab 


for Tub ^ 25 MeV and as 


<7nn = ex P [6.51 exp (T|ab/134) 0 7 ] mb 


( 10 ) 


for Tiab < 25 MeV. 

Note that the energy dependence of the strong interaction cross section is totally contained 
in equation (9). Because of the exponential factor in (5) this energy dependence is rather 
weak as one would expect. 

ITT, electr omagnetic THEORY 

The EM theory has already been discussed extensively (Bertulani and Baur 
1988; Norbury 1989, 1990a,b) and only a few relevant details will be given here. The total 
nucleus-nucleus EM cross section is written as 

a = <j£i + o E2 =| [Nei(E)oei(E) + Ne2(E)c E 2(E)] dE (11) 

where N£i(E) is the virtual photon spectrum (of energy E) of a particular multipolarity i due 
to the projectile nucleus and o E i(E) + a E2 ( E) is the photonuclear reaction cross section of 
the target nucleus. (In principle the above equation should include other EM multipoles, but 
their effect is much less important) A less exact expression is given by WW theory as 


Cww= I Nww(E) [Oei(H) + <Je2(E)] dE 


( 12 ) 


where Nww(E) is the WW virtual photon spectrum. Bertulani and Baur (1988) have shown 
that 


Nww(E) = N E 1 (E) = 12 a X [4 Ro K, - 1 5 2 P 2 <K? - Kj )] (13) 

and 

N E 2(E) = ^r 2 Z 2 a-L [ 2(1-(3 2 )K? + %{2- pVlCoK, - UV(K? - Kg)] (14) 

E ft ^4 l 


with 

% = Ebrnin ( 15 ) 

Y P (he) 

where all of the modified Bessel functions K are functions of In the above equations E is 
the virtual photon energy, Z is the nuclear charge, a is the EM fine structure constant, and 
b m in is the minimum impact parameter, below which the collision occurs via the Strong 
interaction. Also (3 = — and y = = where c is the speed of light and v is the speed of 

C V 1- p 2 

the cosmic ray. The minimum impact parameter is given by 


= b c + 


ft ap 

2y 


(16) 


where 


(17) 


ao = 


Zp Zj e 2 
mo v 2 


allows for deviation of the trajectory from a straight line (Aleixo and Bertulani 1989). 


In equation (1 1) the photonuclear cross sections satisfy the following sum rules 


(Bertulani and Baur 1988): 


c E i(E) dE = 60 MeV mb 

A 


(18) 


and 


ja E2 (E)| = F0.22ZA^ ^ (19) 

where F is the fractional exhaustion of this energy-weighted sum rule. The latter expression 
is the sum rule for the isoscalar E2 giant resonance. The isovector E2 resonance is 
ignored as it decays mainly by 2-nucleon emission (Bertulani and Baur 1988). 

IV. ELECTROMAGNETIC PARAMETRIZAT IQM 

Because the photonuclear cross sections Oei(E) and cte 2 (E) are Lorentzian 
shaped, they behave somewhat like delta functions. The integrals of equation (1 1) can be 
approximated by taking N E i(E) and N E2 (E) outside the integrals as (Bertulani and Baur 


1988): 


a * Nei(Egdr) oei(E) dE + Ne2(Egqr) Eqqr <Je2(E)^| 


( 20 ) 


and the integrals are evaluated using the sum rules in equations (18) and (19). In the above 
equation Eqdr and Egqr are the central energies of the El and E2 photonuclear cross 
sections given by (Westfall, Wilson, Lindstrom, Crawford, Greiner and Heckman 1979) 

Ecdk - he [ (it U - £) ]' n (21) 

with 

u=^f A' 1/3 (22) 

and 

Ro = r 0 A 1/3 (23) 

where e = 0.0768, Q' = 17 MeV, J = 36.8 MeV, r 0 = 1.18 fm, and m* is 7/10 of the 
nucleon mass. Note that other expressions for Egdr such as 80 A 1/3 (Bertulani and Baur 
1988) provide very inaccurate results for light nuclei. Equation (21) is accurate for all mass 
regions. The central energy of the E2 resonance is simply 

Eg< 3 R=-^- MeV (24) 

A l/3 

In addition the fractional exhaustion of the Energy-Weighted Sum Rule in equation (19) is 
given by (Bertrand 1976) 


f = 0.9 for A > 100 


= 0.6 for 40 < A < 100 

= 0.3 for 40 < A (25) 

Finally, to obtain the reaction cross section for proton or neutron removal the 
above cross sections must be multiplied by the proton or neutron branching ratios. The 
proton branching ratio has been parameterized by Westfall, Wilson, Lindstrom, Crawford, 
Greiner and Heckman (1979) as 

g p = min [ f , 1 .95 exp( - 0.075 Z) ] (26) 

where Z is the number of protons and the minimum value of the two quantities in square 
brackets is to be taken. Assuming that only single nucleon removal occurs, the neutron 
braching ratio is 

gn = 1 - g P ( 27 ) 

For light nuclei however the following branching ratios are used instead of equation (26) 


gp - 0.5 

for 

Z < 6 


= 0.6 

for 

6 < Z < 8 


= 0.7 

for 

8<Z< 14 

(28) 


Lastly, an adjustable parameter (the only one in the whole parametrization!) is introduced 
as Xd = 0.25 where 


bmin = (1 + *d) b c + 


(29) 


in place of equation (16). 


Finally, if one is interested in a very quick calculation for estimation purposes we shall 
write down an approximate "pocket" formula which does not require the evaluation of the 
Bessel functions in (13) and (14). Using the low and high frequency approximations for 
the dipole photon spectrum (Jackson, 1975) and ignoring quadrupole effects, equation (20) 
can be written approximately as 


= -J— ^Z 2 aJ-[ln (•5~~T~^ ) ' i P 2 ] for Egdr < 


yhc(3 


Egdr 7t 


r 


Egdr b, 


« Z 2 a -L ( 1 - L p 2 )exp(-2EGDRb mm /7hcp) for E G dr ^ 


Egdr 




yhc(3 


■'min 


(30) 


This EM formula, combined with the Strong interaction parameterization, gives a very 
simple "pocket" formula which may also be useful in complicated versions of transport 
codes that have CPU time at a premium. However to get a good fit to data one must use 


x<i= -0.1 in (30). 


V. RESULTS AND CONCLUSIONS 


The cross section parametrizations are compared with the existing nucleus- 
nucleus experimental data in Tables I - III. It can be seen that the overall agreement is 
extremely good for a very wide variety of projectiles, targets and energies. There are 
however a few notable discrepancies particularly for ^Au targets in Table II. It should be 



noted however that these discrepancies are not due to the parametrization per se. Similar 
discrepancies are observed in comparisons between the original theory and experiment 
(Norbury 1989, 1990a, 1990b, Norbury and Townsend 1990, Benesh, Cook and Vary 
1989, Hill, Wohn, Schwellenbach and Smith 1991). It is not clear whether these 
discrepancies are due to theoretical or experimental problems and their resolution is a matter 
of ongoing research. 

In summary a parametrization of single nucleon removal cross sections for 
nucleus-nucleus collisions has been developed which accurately reproduces the 
experimental data for a wide range of nuclear species and energies. Future work will be 
devoted to few nucleon removal. Combining this with the many nucleon removal 
parametrizations (Wilson, Townsend and Badavi 1987b) and the proton-nucleus 
parametrizations (Letaw, Silberberg and Tsao 1983; Silberberg and Tsao 1973, 1990) 
provides a very useful parametrization of arbitrary cosmic ray species interacting with an 
arbitrary medium. 
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Table I Electromagnetic (EM) Cross Sections for single neutron and single proton 
removal. ofj£ t are the experimental EM cross sections from Olson et al 1981, Heckman 
and Lindstrom 1976, Barrette et al 1990 and Hill 1988. c|£?am is the parameterized EM 
cross section discussed in the text. Values in parentheses use the EM pocket formula. 


Projectile 

Target 

Tiab 

Final 

oEM 

°expt 

t^param 



(GeV/N) 

State 

(mb) 

(mb) 

12c 

Pb 

2.1 

iic 

51 ±18 

48 (57) 

12 C 

Pb 

2.1 

Hb 

50 ±25 

72 (86) 

12 c 

Pb 

1.05 

11c 

39 ±24 

26(17) 

12C 

Pb 

1.05 

11b 

50 ±25 

39 (26) 

16 0 

Pb 

2.1 

15 0 

50 ±24 

71 (85) 

16 0 

Pb 

2.1 

15 n 

96 ±26 

106(127) 

12 C 

Ag 

2.1 

He 

21 ± 10 

20 (24) 

12 c 

Ag 

2.1 

Ub 

18 ± 13 

30 (36) 

12C 

Ag 

1.05 

lie 

21 ± 10 

12(12) 

12 c 

Ag 

1.05 

11b 

25 ± 19 

18(18) 

16 0 

Ag 

2.1 

15 0 

26 ± 13 

29 (35) 

16 0 

Ag 

2.1 

15 n 

30 ± 16 

43 (53) 

12C 

Cu 

2.1 

He 

10 ±7 

9(11) 

12 c 

Cu 

2.1 

11b 

4 ± 8 

13(16) 

12 c 

Cu 

1.05 

He 

9 ± 8 

6(6) 

12 c 

Cu 

1.05 

11b 

5 ± 8 

8(10) 

16 0 

Cu 

2.1 

15 0 

9 ± 8 

13(16) 

!6 0 

Cu 

2.1 

15 n 

15 ± 8 

19 (23) 


Table I continued 


Projectile 

Target 

Tiab 

(GeV/N) 

12 C 

A1 

2.1 

12 C 

A1 

2.1 

12 C 

A1 

1.05 

12 C 

A1 

1.05 

16 o 

A1 

2.1 

16 0 

A1 

2.1 

12c 

C 

2.1 

12 C 

C 

2.1 

12 C 

C 

1.05 

12C 

C 

1.05 

16 0 

C 

2.1 

16 0 

C 

2.1 

18 0 

T i 

1.7 

18 0 

Ti 

1.7 

O 

oo 

Pb 

1.7 

oo 

o 

Pb 

1.7 

18 0 

U 

1.7 

o 

00 

U 

1.7 

32 S 

197 Au 

200 

28 Si 

27 

A1 

13.7 

28 Si 

27 

A1 

13.7 

28 Si 


13.7 

28 Si 

* 2 °Sn 

13.7 

28 Si 

208 pb 

13.7 

28 Si 

2 08pb 

13.7 


Final 

FM 

u expt 

^pwam 

State 

(mb) 

(mb) 

11 C 

0 ± 5 

2(3) 

1 l B 

0 ± 5 

3(4) 

ll C 

1 ± 6 

2(2) 

“B 

1 ±7 

2(3) 

15 0 

0 ± 5 

3(4) 

15 n 

-1 ± 9 

5(6) 

lie 

-2 ±5 

KD 

n B 

-1 ± 4 

KD 

ll C 

-2± 5 

0(1) 

ll B 

-2 ±5 

KD 

15 0 

-1 ± 4 

KD 

15 n 

-1 ± 4 

KD 

l7 0 

8.7 ± 2.7 

8(10) 

17 n 

-0.5 ± 1.0 

12(15) 

* 7 0 

136 ±2.9 

69 (79) 

l7 N 

20.2 ± 1.8 

103(118) 

1 7 0 

140.8 ± 4.1 

82 (92) 

n N 

25.1 ± 1.6 

123 (138) 

196 Au 

1 120 ± 160 

1274(1297) 

IP 

37 ±5 

25 (28) 

In 

15 ±4 

12(13) 

IP 

313 ± 4 

325 (370) 

In 

136 ±6 

151 (172) 

IP 

743 ± 27 

822 (942) 

In 

347 ± 18 

383 (438) 


Table II Total ( = EM + Nuclear) Cross Sections for single neutron removal. aj°p t are the 
experimental total cross sections from Hill, Wohn, Winger and Smith 1988, Hill, Wohn, 
Winger, Khayat, Leininger and Smith 1988, Hill, Wohn, Schwellenbach and Smith 1991, 
Smith et al 1988 and Loveland et al 1988. OpS-’am is the parameterized total cross section 
discussed in the text Values in parentheses use the EM pocket formula. 


Projectile 

Target 

T(ab 

Final 

—Tot 

°expt 

fcTot 

'-'par am 



(GeV/N) 

State 

(mb) 

(mb) 

l2 C 

238u 

2.1 

237y 

173 ±22 

191 (195) 

20 Ne 

238u 

2.1 

237y 

192 ± 16 

286 (300) 

12 C 

^97 Au 

2.1 

196 Au 

178 ±7 

172 (175) 

20 Ne 

197 Au 

2.1 

196 Au 

268 ± 11 

249 (260) 


197 Au 

1.8 

196 Au 

463 ± 30 

458 (491) 

56 p e 

197 Au 

1.7 

^96 Au 

707 ± 52 

748 (812) 


197 Au 

1.26 

^96 Au 

2130 ± 120 

2187 (2295) 

139u 

197 Au 

0.15 

196 Au 

765 ± 48 

729 (883) 

238u 

197 Au 

0.96 

196 Au 

3440 ±210 

3997 (3486) 

16 0 

197 Au 

60 

^96 Au 

400 ±20 

383 (389) 

16q 

197 Au 

200 

^96 Au 

560 ±30 

458 (462) 

12 C 

89y 

2.1 

88y 

1 15 ± 6 

117(119) 

20 Ne 

89y 

2.1 

88y 

160 ±7 

148 (154) 

40at 

89y 

1.8 

88y 

283 ±11 

223 (240) 

56 Fe 

89y 

1.7 

88y 

353 ± 14 

319(351) 

12 C 

59 C o 

2.1 

58 C o 

89 ±5 

99(101) 

20 Ne 

59 C o 

2.1 

58 C o 

132 ±7 

119(122) 

56 Fe 

59 Co 

1.7 

58 C o 

194 ±9 

212 (229) 

139La 

59 C o 

1.26 

58 Co 

450 ± 30 

433 (461) 


Table II continued 


Projectile 

Target 

Tlab 



(GeV/N) 

l2 C 

12C 

2.1 

20 Ne 

12C 

1.05 

56 Fe 

12 C 

1.7 

139La 

12 C 

1.26 

28 S i 

12C 

13.7 


Final 

—Tot 

°expt 

rrTot 
u par am 

State 

(mb) 

(mb) 

He 

61 ± 1 

65 (66) 

He 

78 ±2 

73 (73) 

He 

94 ±2 

100(101) 

He 

148 ±2 

134(135) 

He 

73.5 ± 3.5 

85 (85) 


Table ITT Nuclear Cross Sections for single neutron and single proton removal. c^“p t are the 
experimental nuclear cross sections from Fig. 4 of Barrette et al 1990. is the 

parameterized nuclear cross section discussed in the text 


Projectile 

Target 

Tlab 

Final 

_ nuc 
°expt 

Oparam 



(GeV/N) 

State 

(mb) 

(mb) 

28 S i 

Al 

13.7 

lp 

140 ± 14 

87 

2»Si 

Al 

13.7 

In 

100 ± 10 

87 

28 S i 

Sn 

13.7 

IP 

220 ± 22 

120 

28 Si 

Sn 

13.7 

In 

145 ± 15 

120 

28 Si 

Pb 

13.7 

IP 

300 ±30 

136 

28 S i 

Pb 

13.7 

In 

180 ± 18 

136 



